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SECTION  I 
INTRODUCTION 

A.  MOTIVATION 

Since  the  advent  of  optical  and  infrared  detection  and  tracking 
devices,  the  problem  of  static  and  dynamic  cloud-free  1 ine-of-sight  ICFLOS) 
has  become  increasingly  important.  In  1975,  Nelson  and  Wetherbe  (Ref.  I) 
did  a  first  approximation  study.  They  concluded  that  "it  would  seem  that  the 
probability  of  having  a  cloud-free  1 ine-of-sight,  either  instantaneous  or  for 
as  short  a  period  as  one  second,  will  be  low  even  in  fair  weather  conditions". 
They  also  went  on  to  remark  the  following: 

"To  solve  a  dynamic  CFLOS  problem  one  requires  a  detailed  knowledge  of  the 
initial  positions  and  subsequent  velocities  of  the  two  points  between  which  the 
1 ine-of-sight  is  to  be  assessed,  detailed  knowledge  of  the  boundary  conditions 
related  to  the  space/time  dimensions  of  the  problem,  and  detailed  knowledge  of 
the  cloud  field  characteristics,  including  the  3-dimensional  geometry  and  spatial 
distribution  of  the  individual  cloud  elements.  The  dynamics  of  a  problem  might 
be  definable  through  a  knowledge  of  the  operational  characteristics  of  the  air¬ 
craft  involved.  The  boundary  conditions  and  time  increments  to  be  used  in  the 
1 ine-of-sight  assessment  might  stem  from  the  characteristics  of  the  electro- 
optical  system." 

Thus,  it  became  apparent  to  Nelson  that  a  second  level  of  approximation  was 
needed.  This  led  to  the  development  of  CLOUDS. 

B.  THE  PRESENT  MODEL 

The  second  level  of  approximation  (or  sophistication)  came  into  the 
problem  with  the  development  of  CLOUDS.  This  program  uses  a  three-dimensional 
cloud  field  consisting  of  rotated  elliptical  cylinders  and  rotated  truncated 
ellipsoids.  The  trajectory  generation  capability  in  CLOUDS  consists  of  the 
standard  direct  intercept  and  chase  models  (Section  II.B.l,  II. B. 2),  but 
virtually  any  scenario  can  be  input  through  the  use  of  the  Trajectory  Analysis 
Program  (TAP).  The  statistical  calculations  performed  in  CLOUDS  have  their 


1.  Nelson,  R.  J.  and  Wetherbe,  M.  B.  Some  Aspects  of  Estimating  the  Probability 
of  Cloud-Free  Lines-of-Siqht  in  Dynamic  Situations ,  XFCTAC  TN  76-2.  U.  S.  Air 
Force  Environmental  Technical  Application  Center,  Scott  Air  Force  Base.IL., 
March  1976. 
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data  base  in  CFLOS  for  a  given  time  between  one  point  on  the  defender's 
trajectory  and  one  point  on  one  of  the  aggressor's  trajectories.  Thus,  CLOUDS 
offers  extensive  scenario  generation  flexibility  through  a  three-dimensional 
solution  to  the  CFLOS  problem. 

This  report  is  designed  to  be  an  in-depth  User's  Manual  for  the  CLOUDS 
program.  Examples  of  input  decks  and  the  resulting  output  are  described 
(Section  III  and  IV).  The  overall  program  structure  is  discussed  in  Section  III 
and  crucial  mathematcal  developments  described  in  Appendices  A,  B,  and  C.  The 
CLOUDS  computer  program  is  written  entirely  in  American  National  Standards 
Institute  (ANSI)  FORTRAN  (Ref.  2).  When  compiled  and  run  on  the  AFWL  CDC- 
CY8ER  176  computer,  CLOUDS  requires  53,188  words  (147,704  Octal)  of  core  for 
instructions  and  storage.  Running  time  varies  greatly  with  the  type  of  problem. 
Typically,  for  two  base  heights  (40  percent  cloud  cover  at  each),  the  running 
time  is  under  100  (144  Octal)  seconds. 


2.  ANSI  FORTRAN,  X3-9,  American  National  Standards  Institute,  New  York,  1966. 
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SECTION  II 

COMPUTATIONAL  TECHNIQUE 

The  CLOUDS  computer  program  was  designed  using  modular  construction. 

Modular  construction  of  a  computer  code  offers  many  advantages.  Among  them 
are:  (1)  it  mak.es  debugging  much  easier,  (2)  it  permits  dropping  of  a 

normally  run  portion  of  the  code  with  only  minor  changes  to  the  program 
structure,  and  (3)  it  is  much  easier  for  the  user  to  become  familiar  with 
the  computer  code. 

The  CLOUDS  computer  program  consists  of  two  major  blocks.  The  first  block, 
designated  as  the  CLOUD  FIELD  GENERATION  BLOCK,  reads  the  input  data,  determines 
the  placement  of  the  clouds  via  the  input  data,  and  constructs  the  cloud  field. 
If  desired,  the  cloud  field  is  plotted.  The  second  block,  designated  as  the 
SCENARIO  GENERATION  AND  CFLOS  CALCULATION  BLOCK,  determines  from  the  input  data 
which  method  is  to  be  used  in  scenario  generation  and  performs  CFLOS  and  DCFLOS 
(CFLOS  over  some  finite  time  interval)  calculations.  At  this  point  the  trajec¬ 
tories  of  two  vehicles  can  be  plotted.  The  final  step  in  this  block,  and  in 
the  program,  is  to  output  all  the  relevant  data. 

The  following  is  a  detailed  description  of  each  of  the  modules  except  the 
input  and  output  modules  which  are  discussed  with  examples  in  Sections  III  and 
IV,  respectively. 

A.  CLOUD  FIELD  GENERATION  BLOCK 
1.  Cloud  Placement 

The  CLOUDS  program  can  generate  two  types  of  clouds:  rotated 
elliptical  cylinders  or  rotated  truncated  ellipsoids  (only  the  upper  half  is 
used).  Thus,  the  base  of  either  type  of  cloud  is  an  ellipse.  The  intersection 
between  the  semimajor  and  semiminor  axis  of  this  ellipse  is  known  as  the  cloud 
seed.  Cloud  placement  refers  to  the  placement  of  a  cloud's  seed  in  the  ground- 
based  observer's  reference  frame.  When  clouds  are  placed,  their  seeds  must 
remain  within  a  hemisphere  centered  at  the  ground-based  observer's  origin. 

This  hemisphere  is  known  as  the  sphere  of  cloud  placement.  The  actual  placement 
of  seeds  may  be  performed  in  one  of  two  ways:  overall  placement  or  placement 
by  base  height.  A  description  of  how  to  activate  either  of  these  methods  is 
given  in  Section  III.B.2.  Here  we  will  concentrate  on  the  computational  pro¬ 
cedure  for  each. 
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a.  Overall  (Total)  Placement 

Using  the  user  inputted  base  heights  (minimum  of  1)  and 
total  cloud  cover,  CLOUDS  "randomly"  picks  one  of  the  base  heights  as  the  Z 
coordinate  for  the  seed.  Next,  the  X  and  Y  coordinates  of  the  seed  are  randomly 
picked  so  that  the  seed  remains  within  the  hemisphere  of  the  cloud  placement. 

After  the  seed  is  placed,  the  cloud  is  constructed  from  the  possible  X  lengths, 

Y  widths,  and  Z  thickness  supplied  by  the  user.  This  procedure  is  repeated 
until  the  user  specified  total  cloud  cover  is  reached. 

b.  Placement  by  Base  Height 

Given  the  base  heights  (minimum  of  1)  and  the  percent  of 
desired  cloud  cover  at  each  base  height,  CLOUDS  starts  placing  seeds  in  the  first 
(bottom)  base  height.  Thus,  the  Z  coordinate  of  the  seed  is  "fixed".  The  X 
and  Y  coordinates  of  the  seed  are  found  in  the  same  manner  as  in  the  total 
placement  method.  After  the  specified  cloud  cover  is  reached  in  this  base  height, 

seeds  are  then  placed  in  the  next  base  height,  and  so  on.  When  each  of  the 

base  heights  have  the  specified  cloud  cover,  the  seed  placement  terminates. 

NOTE:  If  one  wishes  to  have  just  one  base  height,  the  two  methods 
of  placement  are  mathematical ly  equivalent.  But  for  efficiency,  the  total 
placement  method  should  be  used. 

2.  Computing  the  Sky  Cover  Between  Two  Base  Heights 

Sky  cover  (the  fraction  of  sky  covered  by  clouds  as  seen  by  a 
ground-based  observer)  is  always  computed  between  an  upper  base  height  (HT2) 
and  a  lower  base  height  (HT1).  This  is  done  to  allow  the  program  to  check  total 
sky  cover  as  well  as  partial  sky  cover  (sky  cover  between  two  base  heights) 
with  the  same  subroutine. 

A  uniform  pattern  of  rays  is  put  up  to  check  the  sky  cover.  These  rays 
emanate  from  the  ground-observer's  origin  and  are  defined  by  each  ray's  direc¬ 
tional  cosines  (Appendix  A).  The  ray  pattern  consists  of  286  rays,  which 
corresponds  to  a  distance  between  rays  of  one-seventh  radian  on  the  unit  hemi¬ 
sphere.  By  changing  the  distance  between  rays  one  can  change  the  number  of  rays 
used.  The  maximum  and  minimum  number  of  rays  that  can  be  used  are  500  and  1, 
respectively.  Statistically,  the  optimum  choice  is  between  150  and  300  rays. 

After  the  ray  pattern  is  set  up,  each  ray  is  compared  with  every  cloud 
until  it  intersects  (is  covered)  a  cloud  or  all  clouds  are  exhausted.  Once 
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a  ray  is  covered  (since  it  is  of  no  furtner  interest),  it  is  "tagged"  and  is 
not  compared  with  any  more  clouds  in  this  and  future  sky  cover  computations, 
when  all  of  the  rays  have  been  compared  for  intersections  with  each  cloud,  the 
sky  cloud  cover  is  computed.  Sky  cover  is  the  ratio  of  the  number  of  covered 
rays  to  the  total  number  of  rays  in  the  pattern. 

3.  Printer  Plots  of  the  Sky  Cover  Between  Two  Heights 

After  the  sky  cover  is  computed,  the  ray  pattern  can  be  trans¬ 
formed  into  a  printer  density  plot.  Points  for  the  plot  are  generated  by  inter¬ 
secting  the  uniform  ray  pattern  with  the  placement  hemisphere.  These  3-space 
(X,  Y,  Z)  points  are  then  transformed  into  2-space  (X,  Y)  points  (Figure  1). 

The  transformation  is  such  that  the  arc  distance  (r)  from  the  pole  to  each 
point  is  preserved  as  a  linear  distance  (r')  from  the  Z-AXIS  to  each  point. 

The  longitude  angle  (o)  for  each  point  is  preserved  as  an  angle  (o')  between 
the  X-AXIS  and  the  projected  distance  r‘. 

Figure  2  is  an  example  of  a  complete  printer  density  plot.  The  periods 
are  rays  that  did  not  intersect  a  cloud  and  the  circles  with  plus  signs  in 
them  are  rays  that  did. 

Like  the  sky  cover  computation,  printer  plots  of  the  sky  cover  are  made 
for  clouds  intersecting  rays  between  two  heights.  In  Figure  2,  the  heights 
are  0,  and  the  maximum  height  a  cloud  top  can  be  away  rom  the  observer's 
origin.  Thus,  the  total  sky  cover  is  displayed  here. 

4.  Computing  the  Ground  Cover 

Ground  cover  is  defined  to  be  the  percent  or  fraction  of  ground 
covered  by  an  orthogonal  projection  of  clouds  (or  parts  of  clouds)  between  two 
heights,  UPHT  and  LWHT,  onto  the  flat  earth  (X,  Y  plane). 

To  compute  the  ground  cover,  the  X.  Y  plane  is  first  set  up  so  that  the 
circle  of  cloud  placement  the  sphere  of  cloud  placement  when  projected  onto 
a  plane  is  a  circle)  is  fully  contained  within  the  grid. 

Check  points  are  the  point'  w^ore  the  X  and  Y  parallel  grid  lines  inter¬ 
sect.  Starting  in  the  upper  left  hand  corner  of  Figure  3,  CLOUDS  works  with 
the  check  points  in  horizontal  rows  until  the  lower  right  corner  is  reached. 
During  processing  CLOUDS  determines  whether  a  given  check  point  is  within  the 
placement  circle.  If  it  is  not,  the  check  point  is  rejected  and  the  program 
moves  on  to  the  next  one.  If  the  check  point  is  within  the  cloud  placement 
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circle,  the  clouds  between  the  two  given  heights  are  inspected  to  see  if  one 
covers  that  check  point.  This  inspection  process  continues  until  a  cloud 
is  found  to  cover  the  check  point,  or  all  the  clouds  are  exhausted.  After 
each  check  point  has  been  examined,  the  ground  cover  can  be  computed.  It  is 
simply  the  number  of  covered  check  points  divided  by  the  number  of  check 
points  within  the  circle  of  cloud  placement. 

5.  Printer  Plots  of  the  Computed  Ground  Cover 

When  computing  the  ground  cover,  CLOUDS  works  across  one 
horizontal  row  at  a  time.  Each  check  point  in  the  row  is  flagged  as  to  whether 
it  is  out  of  the  circle  of  cloud  placement,  in  the  circle  but  clear,  or  lastly, 
in  the  circle  and  covered  with  a  cloud.  The  out-of-range  check  points  are 
transformed  into  blanks;  the  in-range  clear  points  are  transformed  into  periods; 
and  the  covered  check  points  points  are  transformed  into  circles  with  plus 
signs  in  them.  When  a  row  is  completed  it  is  written,  if  desired,  into  the 
output  file.  The  points  in  the  next  row  are  then  transformed  and  so  on. 

Figure  3  is  a  complete  plot  of  the  same  cloud  field  as  in  Figure  2. 

8oth  are  total  sky  views,  that  is,  all  the  clouds  in  the  field  are  repre¬ 
sented.  Notice  (by  comparing  Figures  2  and  3)  that  even  with  the  same  cloud 
field  the  total  sky  cover  and  total  ground  cover  are  not  the  same. 

6.  META  Plot  of  the  Cloud  Field 

At  AFWL  there  exists  a  system  dependent  plot  package  known  as 
METAPLOT. *  CLOUDS  utilizes  this  facility  to  give  the  user  a  three-dimensional 
representation  of  the  cloud  field  and  trajectories  (Figures  4  and  5). 

The  input  data  that  CLOUDS  needs  to  produce  the  META  Plot  is  the 
(X,.  Yj,  Z,)  position  of  the  viewer's  eye  and  the  (X.,,  Y2 ,  Z^.)  position  of 
the  center  of  the  object  to  be  viewed  (Figure  6).  Both  of  these  points  are 
referenced  from  the  ground-based  reference  frame.  Default  values  are  center 
=  (0,  0,  19000),  viewer  =  (90000,  75000,  90000).  The  final  input  parameter 
CL0U0S  needs  to  produce  the  META  Plot  is  NPTS,  the  number  of  vertical  lines 
used  in  drawing  the  clouds.  The  default  value  for  NPTS  is  twelve.  The 
METAPLOT  output  is  placed  on  TAPE99 

NOTE:  If  CLOUDS  is  executed  on  any  computer  other  than  those  at 

AFWL,  the  META  Plot  of  the  cloud  field  should  not  be  activated. 


^Additional  information  on  METAPLOT  may  be  obtained  by  contacting  Bob  Conley  at 
AFWL/DAP,  Kirtland  Air  Force  Base,  NM,  (505)  264-1307,  Autovon  964-1307 
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B.  SCENARIO  GENERATION  AND  CFLOS  CALCULATION  BLOCK 

1.  Chase  Model 

Given  two  vehicles  (i.e.,  their  initial  conditions)  in  the 
ground-based  observer's  reference  frame,  the  chase  model  solves  the  hawk-pigeon 
pursuit  differential  equation  problem.  This  problem  arises  by  requiring  the 
hawk  (or  aggressor,  A)  always  to  fly  directly  towards  the  pigeon  (or  defender, 
B)  at  all  times. 

Additional  restrictions  on  the  trajectories  are:  (1)  the  initial  con¬ 
ditions  on  B  (i.e.,  his  speed,  azimuth,  and  height)  cannot  be  changed.  Thus,  B 
is  restricted  to  a  straight  line  of  flight.  (2)  A's  speed  cannot  change. 

The  algorithm  used  to  generate  the  chase  model  scenario  data  is  discussed  in 
detail  in  Appendix  C.  It  should  be  noted  here  that  if  A's  speed  is  not 
sufficiently  larger  than  B's,  the  two  vehicles  will  never  meet. 

2.  Direct  Intercept  Model 

The  basic  idea  for  this  type  of  scenario  generation  is  that 
both  the  defender  and  aggressor  fly  in  straight  lines  and  that  the  two  vehicles 
must  intersect  at  some  instant  of  time.  The  additional  restrictions  on  this 
model  are  identical  to  the  restrictions  on  the  chase  model. 

For  an  in-depth  look  at  this  problem  see  Appendix  A  of  Reference  1. 

3.  Scenario  Data  from  TAP 

TAP  is  a  very  sophisticated  scenario  generating  program.  It 
is  beyond  the  scope  of  this  report  to  discuss  TAP  in  detail.  More  information 
about  TAP  is  available  in  Reference  3. 

The  aspect  of  TAP  that  is  within  the  scope  of  this  report  is  its  output. 
TAP  has  two  types  of  output,  one  known  as  TACED  (English  units)  and  the  other 
known  as  TAPED  (SI  units).  The  format  in  both  types  of  output  is  identical 
and  CLOUDS  can  utilize  both. 

When  the  TAP  output  is  used  for  scenario  data,  it  must  have  only  one 
defender  and  may  have  up  to  four  aggressors.  Also,  the  data  must  be  placed 
on  a  file  (TAPE7)  so  *  t  CLOUDS  can  read  it. 

3.  Grant,  J.  N. ,  and  Olds,  R.A. ,  "Trajectory  Analysis  Program  (TAP)" , 

AFWL-TR-77-115,  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base, 

NM,  1977.  (Also  contact  AFWL/LEAPS ,  (505)  264-2995) 
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The  following  is  an  example  of  how  to  prepare  the  TAP  output  (on  a 
Permanent  file  called  TACED4)  for  the  CLOUDS  program  on  a  CDC  CYBER  176 
computer. 

JOB 

ACCOUNT 


ATTACH  (TAPE7 ,  TACED4,  ID  =  DYVMMKS) 
REWIND  (TAPE7) 

LGO. 


In  the  ATTACH  statement,  the  TACED4  permanent  file  is  given  the  logical 
file  name  TAPE7.  The  REWIND  statement  just  assures  that  the  pointers  are 
positioned  at  the  beginning  of  the  file.  The  ATTACH  must  be  positioned  in  the 
input  stream  before  the  LGO  so  that  TAPE7  is  set  up  before  CLOUDS  executes. 

4.  CFLOS  Calculation 

For  CFLOS  calculations  only  one  aggressor  at  a  time  is  con¬ 
sidered.  All  of  the  CFLOS  and  DCFLOS  calculations  are  performed  on  a  given 
aggressor  and  then  if  there  is  more  than  one  the  procedure  is  repeated  for 
each.  The  trajectories  through  the  cloud  field  are  contained  in  two  arrays 
(one  for  the  defender  and  one  for  the  aggressor)  as  discrete  points.  The 
time  between  consecutive  positions  is  .At.  CLOUDS  starts  with  the  first  two 
data  points,  one  on  each  trajectory,  at  time  t  =  0.  The  distance  between 
two  points  is  then  checked  to  see  if  it  is  within  the  user-specified  maximum 
and  minimum  allowed  (Section  III. A).  If  not,  this  data  is  rejected,  and  the 
next  point  is  then  inspected.  If  the  distance  is  within  tolerance,  the  line 
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segment  between  the  two  data  points,  known  as  the  1 ine-of-sight  (LOS)  is 
found.  This  LOS  is  then  compared  with  each  cloud  until  an  intersection 
is  found  or  until  all  of  the  clouds  are  exhausted.  The  latter  is  known 
as  cloud-free  1 ine-of-sight  (CFLOS).  The  next  two  data  points  are  inspected, 
and  so  on,  until  all  of  the  trajectory  points  are  exhausted. 

The  total  percent  of  CFLOS  is  the  number  of  cloud-free  1 ine-of-sights 
(within  range)  divided  by  the  number  of  possible  cloud-free  1 ine-of-sights 
(within  range)  times  100. 

5.  DCFLOS  Calculations 

0CFL0S  is  a  consecutive  CFLOS  over  some  time  interval 
T  =  NAt  where  N  >  2. 

In  computing  a  CFLOS,  CLOUDS  sets  up  a  flag  array  with  one  flag  for 
each  pair  of  positions.  The  flags  are  -  1,  0,  and  1,  and  they  correspond 
to  out-of-range,  in-range  without  a  CFLOS,  and  in-range  with  a  CFLOS, 
respectively. 

To  compute  DCFLOS  for  some  given  times,  T,  which  corresponds  to  NAt, 
or  N  positions  in  the  coordinate  arrays  of  the  vehicles,  the  flag  array  is 
inspected.  The  inspection  starts  at  the  first  nonnegative  array  element 
and  looks  for  continuous  groups  of  Nl's.  In  Figure  7,  for  example.  At  =  0.2  s 
and  we  want  a  DCFLOS  for  0.8  s.  This  means  N  =  4,  so  we  look  for  groups  of 
4  continuous  l's.  There  are  4  of  these  groups  out  of  a  total  of  13  pos¬ 
sible  groups.  Thus,  for  this  example,  the  percent  of  0.8  s  DCFLOS  is 
{(4/13)100  =  30.77  percent!. 

The  times,  T,  for  DCFLOS  calculations  are  computed  by  starting  T  equal 
to  2At  and  incrementing  by  2At.  Thus,  T  takes  on  the  values  of  2At,  4At, 

6At,  etc. 


Figure  7.  LOS  Flag  Array  Sample.  The  -Is  indicate  the  LOS  was  out  of  range, the 
Is  Indicate  the  LOS  was  in  range  and  the  sky  was  cloud  free  (CFLOS), 
and  the  0s  Indicate  that  the  LOS  was  in  range  but  covered  with  a 
cloud . 
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If  the  OCFtOS  for  T  goes  to  zero,  computations  are  stopped.  This  Is 
done  since  any  larger  T  will  also  result  In  a  zero  DCFLOS.  In  Figure  7  any 
times  larger  than  1.4  seconds  will  result  In  a  0CFL0S  of  zero. 

6.  Rotations 

Once  a  cloud  field  has  been  generated  It  would  be  a  tremendous 
waste  to  run  the  scenarios  through  It  just  once.  Also,  If  the  scenario  was  by 
chance  oriented  In  the  proper  way,  the  LOS  could  pass  through  a  group  of  clouds, 
or  conversely,  through  a  corridor  between  the  clouds.  Neither  of  these  results 
would  be  representative  of  the  real  physical  situation.  Thus,  In  order  to 
better  utilize  a  generated  cloud  field  and  to  minimize  the  odd  situations 
the  scenarios  are  rotated  NOROT  number  of  times  about  the  Z-axIs.  The  rotations 
are  set  up  so  that  N0R0T  +  1  rotations  bring  each  scenario  back  to  Its  original 
(unrotated)  configuration. 

7.  Plot  of  the  Trajectory 

If  a  METAPLOT  plot  Is  generated  for  the  cloud  field,  the  initial 
(unrotated)  trajectory  of  the  defender  and  each  aggressor  Is  plotted  In  the 
same  plot  frame.  For  details  about  plotting  with  the  METAPLOT  system,  see 
Section  II. A. 6. 
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SECTION  III 
INPUT  DATA  STRUCTURE 

The  complete  input  data  set  for  the  CLOUDS  computer  program  is  divided 
into  problem  sets.  There  can  be  several  problem  sets  in  a  data  set. 

Each  problem  set  is  a  complete  problem  in  and  of  itself.  It  may  or  may  not  be 
related  to  any  other  problem  in  the  data  set. 

Each  problem  set  consists  of  eight  types  of  Hollerith  computer  cards. 

Each  card  type  has  a  specific  field  structure  that  is  unique  to  that  type 
of  card.  The  fields  on  a  card  are  read  in  under  an  NA1  format  where  N  is  the 
length  of  the  field.  Next,  these  Hollerith  fields  are  translated  into 
numbers.  If  the  field  is  blank  or  contains  alphanumeric  (nonnumeric)  symbols, 
then  it  is  translated  as  a  zero.  This  method  allows  the  user  some  flexi¬ 
bility  in  placing  numbers  within  the  field  (e.g.  ,  all  integers  do  not  have 
to  be  right  justified),  and  also  allows  the  use  of  code  words  instead  of 
numbers. 


The  data  fields  used  for  each  card  type  are  described  in  Table  1. 
(Note:  The  International  System  of  Units  (SI)  are  used  in  the  input  data 
as  well  as  throughout  the  program). 


TABLE  1. 

DESCRIPTION  OF 

HOLLERITH  COMPUTER  CARDS 

Card 

Type 

Card 

Columns 

Variable 

Name 

Description 

1 

1-80 

NHEAD 

A  problem  header  statement.  This 
statement  is  printed  out  on  the  first 
page  of  output. 

2 

1-2 

NU1 

The  number  of  base  heights  1  <  NU1  <  15. 

4-13 

PERCNT 

The  total  percent  of  sky  cover  that  is 
desired,  used  only  for  method  1. 

14-17 

PERCLD 

The  percent  of  cylindrical  clouds  to  be 
placed.  The  percent  of  ellipsoidal 
clouds  to  be  placed  is  then  100%- 
PERCLD. 

18-20 

NU2 

The  number  of  possible  different  cloud 
dimensions  1  <  NU2  <  100. 
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TABLE  1.  CONTINUED 


Card 

li£* 

Card 

Co ) umns 

Variable 

Name 

Description 

3 

1-10 

BASHT(I) 

The  Ith  height  at  which  clouds  are  to  be 
placed. 

11-13 

PERCHT(I) 

The  percent  of  sky  cover  desired  at 
this  base  height,  used  only  for  method  2 

21-30 

OENSTV(I) 

The  maximum  density  of  a  cloud  at  this 
base  height  (not  used). 

31-36 

FUN(I) 

The  name  of  a  function  that  describes 
the  density  in  a  cloud  as  a  function  of 
distance  from  the  center  of  the  cloud 
(not  used). 

NOTE: 

There 

must  be  a  cluster 

with  NU1  cards  of  this  card  type  in  it. 

4 

1-10 

L(I) 

The  X-length  of  the  Ith  cloud  type 
in  the  cloud's  reference  frame. 

11-20 

W(I) 

The  Y-width. 

21-30 

T(l) 

The  Z-thickness. 

31-40 

PHI ( I ) 

The  angle  of  rotation  between  the 
clouds  reference  frame  and  the 
ground  based  observer's  refererence 
frame. 

41-43 

PERCDM(I) 

The  percent  of  clouds  that  should 
have  these  dimensions.  PERCDM(l) 

♦  PERCDM(2) 

PERCDM(NU2)  =  100.  This  sum  of 

100  is  3  must! 

NOTE: 

There 

must  be  a  cluster 

with  NU2  cards  of  this  type  in  it. 

5 

1-3 

NPTS 

The  number  of  vertical  lines  used 
to  draw  the  clouds.  If  this  in¬ 
put  is  not  divisible  by  4,  it  will 
be  made  so.  (Default  is  12.) 

5-14 

VIEWER(l) 

The  X-Coordinate  of  the  viewers 
position  in  terms  of  the  ground 
based  observer's  reference  frame. 

15-24 

VIEWERC) 

The  Y-Coordinate 

25-34 

V I  EWE  R ( 3  3 

The  Z-Coordinate.  Default  for  the 

viewers  position  is 
(90000,  75000.  90000) 
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TABLE  1.  CONTINUED 


Card  Card 

Type  Columns 

Variable 

Name 

Description 

41-50 

CENTER(l) 

The  X-Coordinate  of  the  center  of  the 
plot  frame  in  the  ground  based 
observer's  reference  frame. 

51-60 

CENTER(2) 

The  Y-Coordi nate. 

61-70 

CENTER(3) 

The  Z-Coordi nate.  Default  for  the 
center's  position  is  (0,0,16000). 

71-75 

DUMMY 

This  is  the  flag  that  turns  on  the 
METAPLOT  package.  If  DUMMY  >  0  then 
the  METAPLOT  package  is  turned  on 
(Default  =  0,  or  off). 

76-80 

DUMMY2 

This  flag  turns  on  the  printer  plots. 

If  DUMMY2  >  0,  then  the  printer  plots 
are  turned  on.  (Default  is  0,  or 
on). 

NOTE:  See  Reference  3  for  a 

complete  definition  of  viewer  and  center. 

6  1-5 

DELT 

The  time  increment  between  position 
calculations. 

6-10 

DELTAT 

The  maximum  time  for  which  a  continuous 
CFL0S  (DFL0S)  is  calculated.  (Default 
is  the  duration  of  the  scenario). 

11-20 

DMAX 

The  maximum  range  between  aggressor  and 
defender  for  which  a  CFL0S  is  calculated 
(Default  is  10000). 

21-30 

DMIN 

The  minimum  range  between  aggressor  and 
defender  for  which  a  CFL0S  is  calculated 
(Default  is  100). 

7  1-5 

NOROT 

The  number  of  times  (a  maximum  of  18) 
the  scenario  is  rotated  in  the  cloud 
field.  (Default  is  9). 

6-10 

WHSCN 

This  variable  is  a  flag  for  clouds. 

If  it's  positive,  zero,  or  negative,  the 
direct  intercept  model,  the  chase  model, 
or  data  from  the  TAP  file,  respectively, 
will  be  used  for  scenario  generation. 
(Default  is  the  chase  model). 

NOTE:  The  rest  of  this  card  can  be  neglected  if  the  direct  intercept  or 
chase  models  are  used. 
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TABLE  1.  CONTINUED 

Card  Card  Variable 

Type  Columns  Name  Description 

11-15  NOCRFT  The  number  of  aggressors  on  the  TAP 

data  file.  This  number  plus  the  defender 
is  the  total  number  of  crafts  on  the 
data  file. 

16-20  IDENT  The  number  which  identifies  the 

defender  (ownship)  on  the  TAP  data 
file. 

21-30  UNITS  The  type  of  measurement  units  on  the 

TAP  data  file.  Possibilities  are 
"ENGLISH"  or  "MKS" .  (Default  is 
MKS). 


NOTE:  If  the  scenario  data  is  coming  off  the  TAP  file,  this  card  can  be 
neglected  (i.e.,  left  blank),  but  still  must  be  present. 

8  1-10  DAB  The  initial  LOS  distance  between  the 

defender  (B)  and  the  aggressor  (A). 

11-15  ELEVA  The  aggressor's  initial  elevation  angle 

(in  degrees)  from  the  defenders  line  of 
flight. 

16-20  ATZMA  The  aggressor's  initial  azimuth  angle 

(in  degrees)  from  the  defenders  line 
of  flight. 

21-25  ATZMB  The  defender's  initial  azimuth  angle 

(in  degrees)  from  the  ground  based 
observer' s  X-AXIS. 

26-30  VELA  The  aggressor's  speed  during  the 

scenario. 

31-35  VELB  The  defender's  speed  during  the 

scenario. 

36-40  HT  The  defender's  height  above  the  ground. 


B.  INPUT  OPTIONS 

1.  The  following  is  a  list  of  options  that  CLOUDS  can  execute: 
a.  Cloud  Placement 
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(1) 

Clear 

Sky 

(2) 

Overa 1 

1  Placement 

(3) 

Placement  by  Base  Height 

b.  Plotting 

(1)  METAPLOT  -  AFWL  Dependent 

(2)  Printer  Plots  -  Machine  Independent 

c.  Scenario  Generation 

(1)  Direct  Intercept  Model 

(2)  Chase  Model 

(3)  TAP  File 

d.  Multiple  Runs 

2.  Turning  on  the  Options 

The  cloud  placement  options  are  entered  on  the  type  2  and  3 
cards.  It'  a  clear  sky  is  desired,  columns  4-13  on  the  type  2  card  should  con¬ 
tain  the  code  words  CLEAR  SKY.  It  overall  placement  is  desired,  the  same 
columns  on  the  type  2  card  should  contain  a  number  between  0  and  100.  In  the 
latter,  CLOUDS  will  ignore  columns  11-13  on  each  type  3  card,  which  corresponds 
to  the  desired  cloud  cover  at  each  base  height  (Figure  8).  The  user  can  only 
specify  a  total  cloud  cover  or  cloud  cover  at  each  base  height.  Default  is 
total  cloud  cover.  Thus,  if  placement  by  base  height  is  desired,  columns  4-13 
on  the  type  2  card  should  be  left  blank  and  columns  11-13  on  each  type  3  card 
should  contain  a  number  between  0  and  100  (Figure  91. 

The  plotting  options  are  entered  in  columns  71-80  on  the  type  5  card. 
Columns  1-70  all  have  defaults  and  can  be  left  blank.  If  a  METAPLOT  plot  of 
the  sky  is  not  wanted,  columns  71-75  should  be  left  blank  or  a  negative  value 
placed  there.  The  METAPLOT  package  is  turned  on  by  placing  a  positive  value 
in  columns  71-75.  The  METAPLOT  package  is  available  only  at  AFWL  and  CLOUDS 
will  abort  with  unsatisfied  externals  if  the  user  tries  to  use  METAPLOT  at  any 
other  installation.  The  printer  plots  can  be  turned  off  only  by  placing  a 
negative  value  in  columns  76-80  of  the  type  5  card.  If  this  field  is  left 
blank  or  a  positive  value  is  placed  in  it,  the  printer  plots  will  be  used. 
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Printer  plots  are  machine  independent  and  can  be  used  on  any  computer  with  an 
ANSI  compiler  (Figures  8  and  9). 

One  of  three  different  scenario  generation  routines  can  be  selected  by 
placing  a  negative  value,  a  zero,  or  a  positive  value  in  columns  6-10  of  the 
type  card  7.  A  negative  value  turns  on  the  TAP  file  reading  routine. 
Accompanying  a  negative  value  in  this  field  should  be  the  appropriate  infor¬ 
mation  in  columns  11-30  on  the  same  card  (see  Table  1).  The  type  8  card  can 
then  be  left  blank  (but  not  deleted),  as  it  will  be  read  and  ignored.  If  the 
TAP  program  is  used,  the  scenario  data  generated  by  it  must  be  entered  in  the 

TAPE7  file,  so  that  CLOUDS  can  be  read  from  it  (Figure  9  gives  an  example  of 

this  for  a  CDC  CYBER  176  computer). 

The  chase  and  direct  intercept  models  can  be  turned  on  by  placing  a  zero 
or  positive  value  in  columns  6-10  on  the  type  7  card,  or  by  leaving  the 
columns  blank;  CLOUOS  will  ignore  the  rest  of  the  type  7  card. 

The  type  8  card  must  contain  the  information  described  in  Table  1  for  that 

type  of  card.  Multiple  runs  can  be  executed  by  simply  having  more  than  one 
problem  set  in  the  input  data  set  (Figure  8). 
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SECTION  IV 

OUTPUT  DATA  STRUCTURE 

A.  EXECUTION  OUTPUT 

Much  output  data  are  generated  during  execution  of  the  cloud  placement 
phase  of  the  CLOUDS  computer  program.  These  data  (messages)  give  the  user 
a  good  idea  of  what  went  on  during  cloud  placement.  Should  a  fatal  error  occur 
during  processing,  the  messages  will  point  to  the  location  of  the  problem.  The 
order  and  the  type  of  output  messages  depend  on  the  method  of  cloud  place¬ 
ment  being  used  (Section  II.A.l.a  and  b).  The  following  is  an  explanation  of 
the  message  for  each  method  of  cloud  placement. 

PROBLEM  TITLE 

This  is  the  first  program  output  after  the  input  data  has  been 
checked  and  accepted. 

CLOUD  COVER  MESSAGE 

This  message  is  written  after  each  cloud  has  been  placed. 

In  Figure  L0,  starting  from  left  to  right,  the  numbers  correspond  to: 

(1)  The  base  height  for  each  cloud  (HT1), 

(2)  The  upper  height  for  each  level's  sky  cover 
computation, 

(3)  The  number  of  each  cloud  as  it  was  placed,  and 

(4)  The  sky  cover  (between  HT1  and  HT2  for  method  2 
and  total  sky  cover  for  method  1)  for  each  cloud 
as  it  was  placed. 

Note  in  Figure  10,  that  the  first  cloud  given  the  number  12  and  16  was  thrown 
out  and  another  given  its  index  because  the  cloud  placed  first  made  the  sky 
cover  between  4000  m  and  35000  m  greater  than  the  desired  tolerance  of  30 
percent.  If  more  than  10  clouds  are  thrown  out,  a  user  warning  (similar  to 
the  one  shown  in  Figure  11)  is  written  out  and  placement  of  clouds  in  the  next 
level  commences  (method  2)  or  the  sky  is  considered  completed  (method  1  or 
method  2  when  clouds  are  being  placed  in  the  uppermost  base  height). 

BASE  HEIGHT  COMPLETE  MESSAGE 

In  method  2,  when  the  cloud  cover  at  a  base  height  has  been  com¬ 
pleted  (all  clouds  at  that  base  height  have  been  placed),  a  message  (similar  to 
the  one  shown  in  Figure  12)  is  written  out. 
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Figure  10.  Cloud  Placeaent  Messages  Detailing  the  Sky  Cover  After  Each  Cloud 
Is  Placed. 
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RESIDUAL  CLOUD  COVER  ME'- SAGE 

Assuming  again  that  method  2  is  being  used,  when  all  clouds  at  a 
base  height  are  placed,  the  cloud  cover  tor  the  next  base  height  is  computed. 
Since  no  clouds  have  been  placed  in  the  new  base  height,  the  cloud  cover  computed 
here  is  from  clouds  coming  up  from  below.  This  number  is  written  to  the  output 
file  with  an  appropriate  message  (Figure  12). 

SKY  COMPLETE  MESSAGE 

When  the  total  sky  cover  is  satisfied  for  method  1  or  the  sky 
covers  for  all  base  heights  have  been  completed  for  method  2,  the  fraction  of 
total  sky  cover  is  computed  and  written  to  the  output  file. 

When  method  1  is  selected  for  cloud  placement,  the  order  of  output  is: 

(1)  Problem  title. 

(2)  Cloud  cover  messages  until  the  sky  is  completed.  If  desired, 
printer  plots  of  the  total  sky  cover  is  written  to  the  output 
file. 

(3)  The  cloud  cover  is  computed  between  each  two  base  heights. 

If  desired  printer  plots  of  the  cloud  cover  between  each 
two  base  heights  are  written  to  the  output  file. 

For  method  2,  the  messages  start  out  quite  similar  to  those  for 
method  1  but  diverge  rapidly.  The  order  of  output  for  method  2  is: 

(1)  Problem  title. 

(2)  Cloud  cover  messages  are  generated  until  a  base  height  is 
completed.  If  printer  plots  are  being  produced,  the  sky 
cover  printer  plot  for  this  base  height  is  written  to  the 
output  file. 

(3)  A  complete  message  is  written  out  for  each  base  height 
as  it  is  completed. 

(4)  If  the  uppermost  base  height  has  not  been  reached,  the 
residual  cloud  cover  is  computed  for  the  next  higher  base 
height  and  printed.  Number  2  messages  are  generated  until 
all  of  the  base  heights  are  completed. 

(5)  The  total  sky  cover  is  computed  and,  if  desired,  the  total 
cloud  cover  printer  plot  is  generated. 

At  this  point  methods  1  and  2  become  similar  again;  both  go  on  to 
the  ground  cover  determinations  (Section  1 1. A. 4  and  51,  which  consist  of: 

(11  total  ground  cover,  and  (2)  ground  cover  due  to  clouds 
between  base  heights. 
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Figure  13.  Scenario,  Rotation  and  CFIOS  Generation  Messages.  Each  Message  1i 
Printed  after  the  Activity  It  Describes  Is  Coapleted. 


Table  2.  Descriptive  Date  for  Each  Cloud  Placed. 
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Table  3.  First  Aggressor  CFLOS  Calculations 
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During  each  ground  cover  computation  a  printer  plot  can  be  generated  and  written 
to  the  output  tile. 

The  cloud  generation  phase  of  execution  has  now  been  completed,  and  the 
scenario  generation  and  CK  LOS  calculations  begin.  I  he  messages  tor  this  phase 
are  self-explanatory  (Figure  1.1). 

B.  POST  EXECUTION  OUTPUT 

This  phase  of  output  displays  the  important  data  generated  in  a 
problem  set.  The  post  execution  output  (Table  2)  contains  all  of  the  descriptive 
data  tor  each  cloud.  Table  0  gives  the  user,  along  with  the  execution  output 
cloud  cover  messages,  useful  information  about  each  cloud  Next,  the  input 
data,  as  translated  and  used  by  CLOUDS,  are  listed  (Figure  1*1).  For  a  complete 
description  of  the  input  data,  see  Section  111.  The  input  data  are  followed 
by  a  table  containing  a  summary  of  the  sky  and  ground  covers,  both  total 
coverage  and  coverage  between  base  heights. 

The  results  of  the  CF LOS  calculations  are  then  written  in  the  form  of 
two  tables  for  each  aggressor  The  first  table  (Table  .))  contains  the  per¬ 
cent  of  total  Cl  LOS  and  the  longest  continuous  CFLOS  in  seconds  for  each 
rotation.  Ihe  third  column  in  this  table  contains  the  number  of  times  the 
longest  continuous  CFLOS  occurred.  The  next  table  (Table  4)  contains  the  per¬ 
cent  of  DCFLOS  for  different  times  and  rotations.  The  times  for  each  DCFLOS 
are  found  by  taking  the  integer  under  the  heading  T1MI  and  multiplying  it  by 
the  time  increment 

After  the  DCFLOS  table  for  the  last  aggressor  is  listed,  the  preset  random 
number  generator  is  printed.  This  number  is  the  time  that  CLOUDS  started 
computing.  If  the  user  desired  to  duplicate  a  run  exactly,  this  number  must  be 
placed  in  the  variable  1 T l ME’  before  the  call  to  RANSE  1  on  the  main  program. 

The  format  for  setting  ITIMF  is: 

ITIME  =  IDHbHH.MM. SS. 

In  the  output  message  for  the  random  number  generator,  the  last  period  is 
part  of  the  preset  variable  ITIME,  not  punctuation,  and  must  be  in  the 
setting  for  IT  IMF. 
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APPENDIX  A 
THE  RAY  PATTERN 

It  will  become  evident  in  Appendix  B  that  a  very  important  tool  is  needed 
to  check,  the  sky  cover  after  each  cloud  is  placed.  This  tool  is  a  "UNIFORMLY" 
distributed  ray  pattern.*  The  generation  of  a  ray  nattern  that  is  perfectly 
uniform  is  a  very  interesting  and  quite  involved  statistical  problem.  Since 
the  CFLOS  problem,  rather  than  the  perfectly  uniform  ray  pattern  problem, 
was  our  main  concern  an  approximation  to  the  perfectly  uniform  ray  pattern 
was  developed. 

The  approximation  is  based  on  the  rather  simple  idea  that  the  intersections 
of  the  rays  and  unit  hemisphere  should  be  equally  spaced  on  equally  spaced 
latitude  lines.  To  implement  this,  the  following  algorithm  was  used: 

Let  0  be  the  arc  distance  between  latitudes.  Then  the  equally  spaced 
latitudes  can  be  generated  by  starting  at  the  zenith  (or  2-axis)  and  incre¬ 
menting  along  the  spheres  surface  by  0  (Figure  A-l).  At  each  latitude  the 
points  of  intersection  are  to  be  equally  spaced  so  that  D  must  be  slightly 
adjusted.  Let  NMAX  be  the  number  of  intersections  on  a  given  latitude  and 
DLON  the  increments  between  the  intersections,  then  NMAX  and  DLON  are  given 
by: 

NMAX  =  IFIX(2n*SIN(N*D)/D  ♦  0.5) 

DLON  =  2n*SIN(N*D) /FLOAT (NMAX) 

Where  N  is  the  number  of  the  latitude  line.  Points  of  intersection  are 
placed  on  a  latitude  line  by  starting  on  the  longitude  line  that  lies  in  the 
X-Z  plane  and  incrementing  around  the  latitude  line  by  increments  of  DLON. 

To  make  the  problem  more  uniform  this  process  is  done  on  every  other  latitude 
line  (the  odd  numbered  ones).  For  the  even  numbered  latitude  lines  the  first 
incremental  distance  from  the  X-Z  longitude  line  is  QLON/2,  but  everything  else 
remains  the  same. 

When  generating  a  ray  pattern  it  is  desirable  to  use  directional  cosines 
for  each  ray.  See  Appendix  B  for  a  definition  of  directional  cosine.  Suppose 
we  have  a  point  of  intersection,  P,  as  labeled  in  Figure  A-l.  Then,  the 


*A  uniformly  distributed  ray  pattern  refers  to  a  ray  pattern  that  when  inter 
sected  with  the  unit  hemisphere,  gives  any  two  areas,  with  the  same  size, 
the  same  probability  of  having  a  ray  intersect  the  surface  within  them. 
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Figure  A-1 .  The  Direction  Cosines  (a,  S,  y)  for  a  Ray  Intersecting 
Hemisphere  at  the  Point  P. 
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directional  cosines  for  the  ray  that  intersects  the  unit  hemisphere  to  produce 
P,  are  A,  B,  and  G.  To  find  these  distances  we  need  y  and  0.  These  are  given 
by: 

y  =  N  *  D 

9  =  (M  *  DION  +  BIAS)/r 

where  M  is  the  number  of  the  present  increment  around  the  latitude  line 
and  BIAS  is  given  by: 

10  for  N  =  2K  +  1 

DL0N/2  for  N  =  2K 

Thus,  it  is  apparent  that 
A  =  r  *  cos0 
B  =  r  *  sin0 
G  =  cosy 

Where 

r  =  siny 

Thus,  with  the  above  algorithm  one  can  find  the  directional  cosines  for  a 
UNIFORMLY  distributed  ray  pattern. 
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APPENDIX  B 
INTERSECTIONS 

When  computing  sky  cover  a  number  of  rays  are  sent  out  (Appendix  A)  and 
each  ray  that  has  not  previously  been  covered  is  checked  against  each  cloud 
until  an  intersection  is  found  (thereby  covering  the  ray).  Clouds  are  repre¬ 
sented  by  rotated  truncated  ellipsoids  (only  the  upper  half  is  used)  and 
rotated  elliptical  cylinders.  Thus,  for  computational  speed  it  is  imperative 
that  the  solution  to  these  equations  of  intersection  be  developed  and  written 
in  a  form  best  suited  for  numerical  evaluation.  What  follows  is  such  a 
development.  First,  the  equation  for  the  ray-rotated  elliptical  cylinder 
intersection  in  three-dimensional  space: 

( Xcos<|»  +  Vs i nq>  -  XI)2  .  (-Xsin<j>  +  Ycos<j>  -  Yl)2  _  ,  (B-l) 

Z1  <  Z  <  Z1  ♦  T 

Where  L,  W,  and  T  are  the  length,  width,  and  thickness,  respectively  of  the 
cloud  in  the  primed  reference  frame  (Figure  B-l);  and  $  is  the  angle  of  rotation 
between  the  primed  (cloud  centered)  reference  frame  and  the  observers  refer¬ 
ence  frame.  The  restriction  on  Z  makes  the  rotated  ellipse  into  a  rotated 
elliptical  cylinder. 

The  equations  of  a  line  in  three-dimensional  space  is  given  by: 

X  =  at  +  PI 

Y  =  pt  +  P2  (B-2) 

Z  =  Yt  ♦  P3 

Where 

a  =  cose 

p  =  cos6 

Y  =  cos£ 

Thus  a,  p,  y  are  the  directional  cosines  of  the  ray  originating  at 
(PI,  P2,  P3). 
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Substituting  Equation  B-2  into  Equation  B-l  yields  a  quadratic  equation 
in  t  as  follows: 

(Ata2  ♦  A2p2  ♦  A3y2  +  A4a|i)  t2  +  (Bi<i  ♦  B2p  +  B3y  ♦  B4aPl 

♦  BsaP2  ♦  B6pPl  B7pP2  ♦  B9yP3)  t  ♦  (C,P12  ♦  C2P22  +  (B-3) 

C3P32  ♦  C4P1P2  ♦  CSP1  ♦  CSP2  +  C7P3  ♦  C9)  =  0 

Where 

kx  =  W2cos2$  ♦  L2sin2$ 

A2  =  W2sin2#  L2cos2$ 

A3  =  2cos<J>sin<Ji(W2  -  L2) 

A4  =  0 

Bj  =  -2(X1W2cos4>  -  YlL2sin$) 

Bo  =  -2(XlW2si n$  +  YlL2cos<t> 


B3  =  0  -  A4 

B4  =  2(W2cos24>  ♦  L2sin2<j>)  =  2Aj 

B5  =  2cos$sin$(W2  -  L2)  =  A3 

B6  =  2cost>sin<>(W2  -  L2)  =  A3 

B7  =  2(W2sin2$  ♦  L2cos2iJ>)  =  2A2 

B8  =  0  =  2A4 

Ct  =  W2cos2<>  ♦  L2 s i n2<j»  =  Ai 

C2  =  W2sin2$  ♦  L2cos2$  =  A2 

C3  =  0  =  A4 

C4  =  2sin$cos»KW2  -  L2)  =  A3 

C5  =  -2(W2Xlcos<>  -  L2YlcosO)  =  Bx 

C6  =  - 2(W2 Xls i n4>  ♦  2L2Ylcos<}>)  =  B2 

C7  =  0  =  B3 


C9  =  L2(X1  -  W)  (Y1  +  W)  ♦  W2X12 

In  developing  the  ray-rotated  truncated  ellipsoid  intersection  a  full 
(Nontruncated)  rotated  ellipsoid  will  be  used.  Thus,  after  intersections  are 
found  they  must  be  checked  to  make  sure  they  lie  in  the  upper  half  of  the  full 
ellipsoid.  The  equation  for  a  rotated  ellipsoid  in  three-dimensional  space  is 
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(Xcoso+Ysinq-Xl)2  +  (-Xsin^Ycos^-Yl)2  +  (Z-Zl)2  _  1  (g.4) 

L2  W2  T2 

Substituting  equation  (B-2)  into  equation  (B-4)  yields  a  quadratic  equation 
in  t  equivalent  to  equation  (B-3),  where  the  A's,  B's  and  C's  are  defined  as 
follows: 


A!  =  T2(W2cos2$  +  l2sin2i)>) 

=  T2(W2sin2<J>  +  l2cos2$) 

A3  =  2T2cos<$sin<(>  (W2  -  L2) 

A4  =  W2L2 

B!  =  -2T2(XlW2cos0  -  YlL2sin<f) 

B2  =  -2T2(XlW2sin$  +  YlL2cos<)>) 

B3  =  -2Z1L2W2 


B4  =  2T2(W2cos2$  +  l_2sin2q)  =  2At 

Bs  =  2T2cos4)sin<ti(W2  -  L2)  =  A3 

Bs  =  2T2cos<J>sin<)>(W2  -  L2)  =  A3 

87  =  2T2(W2sin24>  +  L2cos2$)  =  2A2 

Bs  =  2W2L2  =  2A4 

Ci  =  T2(W2cos2q  +  L2sin2(j>)  = 

C2  =  T2(W2sin2t)>  +  L2cos2«f> )  =  A2 

C3  =  W2L2  =  A4 

C4  =  2T2cos4»sin4>(W2  -  L2)  =  A3 

C5  =  -2T2(XlW2cosq  -  YlL2sin<(>)  =  B1 

C6  =  -2T2(XlW2sin<J>  ♦  Yll2cos4>  =  B2 

C7  =  ~2Z1W2L2  =  B3 


C8  =  T2W2  (XI  -  L)  (XI  +  L)  +  t2(T2Yl2  +  W2Z12) 
Now  let  equation  (B-3)  be  represented  by: 


(B-5) 


(B-7) 


At2  ♦  Bt  ♦  C  =  0  (B-8) 

This  quadratic  equation  yields  itself  very  quickly  to  analysis. 

If: 


B2  -  4AC  <0 
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then  the  roots  of  equation  (B-8)  are  complex  (^0)  or  double  (=  0)  and  both  are 
counted  as  nonintersections.  If  equation  (B-9)  does  not  hold,  then  the  inter¬ 
sections  are  given  by: 


It  and  I2  can  now  be  analyzed  to  determine  if  they  lie  in  the  "cloud" 
region  of  the  rotated  ellipsoidal  cylinder  or  the  rotated  truncated  ellipsoid. 
For  the  rotated  ellipsoidal  cylinder  the  points  of  intersection  must  satisfy: 

(Z1  <  It  <  Z1  +  T)  and  (Z1  <  I2  <  Zl  ♦  T) 
or 

( I !  <  Zl)  and  (I2  >  Zl  +  T) 

For  the  rotated  truncated  ellipsoid  the  points  of  intersection  must  satsify: 
(It  >  Zl)  and  (I2  >  Zl) 
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APPENDIX  C 
SCENARIOS 

The  CLOUDS  computer  program  has  two  models  to  internally  generate  the 
scenarios.  These  models  are  known  as  the  direct  intercept  model  and  the  chase 

model.  A  development  of  the  direct  intercept  model  can  be  found  on  page  16  of 

Reference  3.  What  follows  here  is  a  predictor-corrector  solution  of  the  chase 
model.  See  Section  II.B.l  for  a  description  of  the  chase  model. 

Given: 

(1)  Z  height  of  B,  HT 

(2)  Azimuth  of  B  with  respect  to  the  X-axis  0g 

(3)  Azimuth  of  A  with  respect  to  B's  line  of  flight, 

(4)  Inclination  of  A  (at  t  =  0)  with  respect  to  the  Z  =  HT  plane,  y 

(5)  The  line  of  sight  distance  between  A  and  B  at  t  =  0,  DAB, 

(Figure  C-l). 

(6)  Speed  of  A  and  B,  VA  and  Vg  respectively 

(7)  Time  increment  between  position  calculations,  DELT. 

The  problem  is  to  find  the  path  of  A  so  that  the  tangent  to  A's  path  is 
always  intersects  8. 

To  solve  this  problem  an  iterative  predictor-corrector  scheme  was  used. 
The  steps  of  the  scheme  are  as  follows: 

(1)  Find  the  initial  positions 


XA  = 

0 

DAB*cosycosa 

YA  = 

0 

DAB*cosysi na 

ZAo  = 

DAB*siny  +  HT 

XB  = 

0 

0.0 

II 

o 

CO 

>- 

0.0 

ZB  = 

HT 

where 
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(2)  Compute  the  directional  cosines  for  the  motion  of  both  crafts, 


using: 

«o*  vK  -  XBo)!  *  K  -  VBo)'  *  K  -  ZBo)^ 

The  directional  cosines  for  A  are: 

DXA  =  Va-(XB0  -  XA0)/p0 

DVA  =  V^(YBo  -  YAq)/po 

DZA  =  V.*(ZB  -  ZA  Vp 

A  V  o  oj  Ko 

and  for  B  the  directional  cosines  are: 


DBX  = 


Ycos6B 


using: 


then 


DYB  = 

Vs 

ineB 

DZB  = 

0.0 

(3)  Move 

the 

defender  (B)  one  time  step 

*Bn+l  = 

XB 

n 

+  DXB*DELT 

YVl  = 

YBn 

+  DYB*DELT 

ZBn+l  = 

ZB 

n 

(4)  Compute  the  new  directional  cosines 

for  A 

and  take  the  average 

-A 

XVl)2  *  K  -  VVl)1  * 

(zv 

■  ZBn+l)2 

DXA 

new 

V( 

XBn  -  XBn+l)/pn 

DYA  = 

new 

V( 

VB„  -  ¥Bn.l)'p„ 

DZA 

new 

v( 

-  “n4/pn 

DXA  = 

(dxa 

*  °xV.)/2-° 

DYA  = 

(dya 

♦  DYA  i/2.0 

new/ 

DZA  = 

(DZA 

+  DZAnowV2-0 

new ) 

;  the  "=" 

sign 

in  the  last  three  equations 

is  a 

fortran  "="  sign  and 

is  read  DXA  "becomes"  .... 

(5)  Move  the  aggressor  one  time  step 
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XArvl  =  XAn  *  DXA*DEET 

YAn>l  =  YAn  *  DYA*DELT 

ZAn+l  =  ZAn  +  DZA*DELT 

(6)  Go  back  to  3  and  repeat  the  process  until  the  crafts  are  within 

0. 5  meter. 

The  accuracy  to  which  the  trajectory  of  A  is  calculated  can  be  increased 
by  iterating  steps  three  through  five  100  extra  times  (and  thus  a  time 
increment  of  0ELT/1Q0  must  be  used)  and  storing  every  100th  position  that  is 
calculated. 
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